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O 1. Introduction 

3 : 

£3 ■ Predictions and structure studies of superheavy elements (SHE) have 

been made, since the establishment of the nuclear shell modelpQ. Espe- 
cially in the last decade, elaborate investigations have been performed of 
rN \ shell correction energy and thereby of a possible location of the superheavy 

island in the nuclear chart. Furthermore, not only about the center of the 
island, but also about stability properties of nearby nuclei have been be- 
ing investigated, which is useful for the extension of the nuclear chart in 
heavy and superheavy elements^. On the other hand, studies of nuclear 
reaction mechanisms have not been developed so much, though so-called 
fusion-hindrance was experimentally found to exist in heavy ion fusions and 
inferred to be due to energy dissipation [3|. That is, there is no reliable the- 
oretical framework which enables us to predict fusion probability of massive 
systems and thereby residue cross sections of SHE properly. Thus, which 
combination of incident ions is most promising and what incident energy 
is an optimum is not yet predicted theoretically. Therefore, the fusion ex- 
periments have been performed, based on systematics of data available so 
far Q. 

Based on the reaction theory of the compound nucleus [5], residue cross 
sections are given as follows, 
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£ T = vrX 2 S(2J+l).P f 4 ion -P/ urv , (1) 

where X is the wave length divided by 2 • it and J the total angular momen- 
tum of the system. Pf us ion and -P surv are fusion and survival probabilities, 
respectively. In the present paper, we discuss several difficult problems in- 
herent in synthesis of the superheavy elements with brief explanations of a 
few progresses of our understanding, as well as attempts of realistic calcu- 
lations. 



2. Difficulties Characteristic in Synthesis of SHE 

In order to obtain the fusion probability, we have to take into account 
possible mechanisms for the fusion-hindrance. Otherwise, calculated proba- 
bilities, and fusion cross sections would be unrealistic, as it is the case that 
one uses a transmission coefficient of an optical model or a barrier penetra- 
tion factor as the fusion probability. As for possible origins of the hindrance, 
two mechanisms are proposed. One is dissipation of incident energy in the 
course of two-body collisions and thus probability for the system to overcome 
the Coulomb barrier is reduced [Jj. The other one is dissipation of energy 
of collective motions of the amalgamated system which has to overcome a 
conditional saddle or a ridge line in order to reach the spherical shape, i.e., 
the compound nucleus Thus, the probability for reaching the spherical 
shape is also reduced. It is natural to consider that both exist. In other 
words, the fusion probability Pf us ion consists of two factors; the sticking 
probability P s ti c k of two incident ions after overcoming the Coulomb barrier 
and the formation probability Pf or m of the spherical shape after overcom- 
ing the conditional saddle point, starting from a pear-shaped configuration 
made by the sticking of the incident ions[S]. 

pJ __ pJ pJ 

r fusion -'stick ' -'form 

Since the existence of the saddle point or the ridge line between a pear- 
shape made by the incident ions and the spherical shape is typical in very 
heavy systems, the latter mechanism would be indispensable for the fusion- 
hindrance observed in massive systems, though the former would also play 
a role (Note that in lighter heavy-ion systems the amalgamated shape is 
usually located inside the ridge line, so the system eventually slides down 
to the spherical shape with probability being equal to 1, once the incident 
ions stick to each other, though fluctuations to be discussed below may 
reduce it only slightly). In either mechanism, we have to describe a passing 
over a barrier under energy dissipation, which is not yet well understood 
theoretically jlUj and thereby there is no useful formula ready for practical 
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applications. This is quite contrast to a similar problem, i.e., to fission under 
dissipation, where a famous Kramers formula^ for decay rate is known to 
well describe a process of a system inside a potential pocket leaking over the 
fission barrier. An essential difference is that in the latter the initial system 
is in the quasi-equilibrium in the pocket, while in the former the initial 
state is given by the condition of two incident ions with a given incident 
cm. energy. 

Recently, the present author and his collaborators have proposed a new 
analytic formula for the probability of passing over a parabolic barrier under 
frictional force. We have applied this formula to the problem of passing- 
over a conditional saddle point, and obtained a simple expression for so- 
called extra-push energy which provides a clear understanding of the fusion 
hindrance |12l I13| . This would be an important contribution to study of 
fusion mechanisms and is briefly recapitulated in section 3., but there still 
remains a difficulty in practice. The parabolic shape is usually a good 
approximation for barrier shapes, but in potential landscapes calculated 
with the liquid drop model (LDM) a pocket inside the saddle is very shallow 
in nuclei corresponding to the superheavy elements, as is easily expected 
from the fissility parameter Xf being close to 1. Therefore, the potential is 
expected to be substantially asymmetric around the saddle, and moreover 
a system once passing over the saddle may return back with an appreciable 
probability due to the fluctuation associated with the friction. Of course, 
the probability for return-back to re-separation is reduced if the system is 
cooled down by neutron emissions and restores the shell correction energy 
which makes the pocket deeper. 

For a quantitative prediction, those features should be taken into account 
properly, which is made by numerically solving a Langevin equation 44 . 
For a dynamical description of shape evolutions, we have to solve trajecto- 
ries in a multi-dimensional space of shape parameterization with a realistic 
LDM potential, examples of which are discussed in section 5. But a time- 
dependent shell correction energy due to evaporation of neutrons is not yet 
taken into account in fusion processes. (Since time for fusion process is 
expected to be rather short, this would not cause a serious inaccuracy, but 
is properly done in the calculation of survival probability, as will be dis- 
cussed in sections 6. and 7.) In the approaching phase of passing over the 
Coulomb barrier under friction, we have to take into account a dissipation 
of the orbital angular momentum as well as that of the kinetic energy of 
the radial motion, where a coupling between them is not described by a 
quadratic potential, as is discussed in section 4. Of course, there are many 
other effects which may play a role in the latter process, say, effects of de- 
formations of incident ions, quantum tunneling effects etc., which are not 
yet fully investigated. 
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As for the survival probability, the statistical theory of decay is well 
established for obtaining a probability for the system to survive against fis- 
sion and charged particle decay. But in practice there are ambiguities in 
the physical parameters, i.e., so-called level-density parameter a and the 
shell damping energy which controls restoration of the shell correction 
energy by cooling. Especially, the latter is crucially important, because the 
restoring shell correction energy gives rise to an additional fission barrier 
effectively which controls survival probability, which is discussed qualita- 
tively in section 6. and quantitatively in section 7. In addition, there 
are Kramers |lll I15j and collective enhancement factors [To] in fission decay 
widths to be taken into account. These are briefly discussed and exam- 
ples of realistic calculations on Ca+ actinide targets are presented which 
are made by employing a new statistical code KEWPIE|17j for the survival 
probabilities, in section 6. 



3. Fusion Hindrance and Extra-Push Energy : Parabolic Barrier 

We study a problem of obtaining a probability for passing over a po- 
tential barrier under a frictional force, which originates from interactions 
of the degree of freedom under investigation with a heat bath, i.e., with 
other degrees of freedom. Therefore, there should be a random force associ- 
ated with the friction in accord with the dissipation-fluctuation theorem. If 
we approximate the barrier with an inverted parabolic shape, the equation 
of motion for a coordinate q and its associate momentum p is written as 
follows, 

d ( q \ _ ( 1/m \ / q \ / 



dt V P J K rnuj 2 -13 ) \ p J + { R J ' (3) 

where m denotes the inertia mass, and uj the curvature of the inverted 
parabola. (3 is a reduced friction, i.e., the friction 7 divided by m, while R 
is its associated random force. The random force is assumed to be Gaussian 
and satisfies the flowing properties, 



(R(t)} = 0, (4) 
(R(t)-R(t')) = 2-7-T-£(t-i / ), 

where () signifies an average over all the possible realizations and the last 
equation given in Eq. (4) with temperature T of the heat bath expresses the 
dissipation- fluctuation theorem. Since the equation is linear, one can write 
down a general solution. 
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With this solution an general expression for a distribution function 
W(q,p;t) at any later time t is calculated, starting with the following defi- 
nition, 



W(q,p; t) = (5(q- q(t)) ■ 5 (p - p(t))) {R} , (5) 

where q(t) and p(t) denote a general solution of Eq. (3) and is given by 
a linear combination of their initial value go an d Po with the coefficients 
including parameters j3 and uj. 0{,r} again denotes the average over all 
the possible realizations of R(t). Using the path integral technique, we can 
perform the averaging and obtain the distribution function of the system 
as a Gaussian distribution around the mean trajectory ((q(t)) , (p(£)))[T2j. 
Then, a probability for passing over the barrier is calculated by integrating 
over the whole p-space and the half g-space, and then by taking the limit of 
time t to the infinity, 



Pform 2 




X 




(6) 



where x denotes (3 divided by 2 • uj. K and B denote the initial kinetic 
energy Pq/2^ and the barrier height measured from the initial potential 
energy \i ■ uj 2 ■ q^/2, respectively. 

In order for the probability to be 1/2, the argument of the error function 
should be equal to zero, which means that the mean trajectory just reaches 
at the top of the barrier overcoming the friction. Then, the necessary critical 
kinetic energy K c is given as 

/ / S 2 

K, 



(Vx 2 + l + x) -B, (7) 



where we can see that in the case of no friction, i.e., of x being equal to 
zero, K c = B, which is trivial. It clearly shows that K c is much larger 
than B under the frictional force. If we estimate the first factor in Eq. (7), 
assuming One-Body Wall-and- Window formula (OBM 18 ) for the friction 
7, we obtain about 10, depending on a reasonable choice of values for the 
inertia mass and the curvature of the potential calculated with LDM. This 
gives a simple formula for the extra-push energy, though we should be careful 
in a comparison with experimental data about effective one-dimensional 
quantities for B, p,, u> and Coulomb barrier heights of entrance channel etc. 

Another interesting formula is obtained, which is very useful for synthesis 
of SHE. Residue cross sections are extremely small in SHE, i.e., we are facing 
with the situation where fusion probability is very small. This suggests in 
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our present formulation that the mean trajectory does not reach the top of 
the barrier, even is far before the top, which means that the argument of 
the error function of Eq. (6) is very large. Then, employing an asymptotic 
expansion of the error function, we can obtain a simple approximate formula 
for the formation probability Pr^ m [T3] . 



Pfo 



1 1 




(8) 



• cxp 



where it is interesting to note that there is a factor very similar to Arrhenius 
factor which is typical in thermal activation processes such as nuclear fis- 
sion, neutron evaporation, thermal electron emission from metal, etc. The 
exact Arrhenius factor is obtained in case of a complete damping of the 
relative motion in the approaching phase, as follows. As will be shown 
in the next section, a distribution of the radial momentum at the contact 
point is approximately expressed by a Gaussian as a results of two-body 
collision processes and thus, the formation probability is obtained by a con- 
volution over initial momentum po, which results again in an error function 
of Eq. (6) with K being replaced with the average value K and with the 
associated variance. In case of completely damping, K is equal to zero, and 
the variance is equal to the temperature. Accordingly, the corresponding 
asymptotic expansion gives the exact Arrhenius factor. 



Pfo 



1 1 



exp 



B 
T 



0) 



Since fusion is inverse to fission in reaction directions, one could call this 
as an "inverse Kramers for inula [T§]" . But we should be careful that in the 
thermal activation processes the factor appears in decay rate or in emission 
rate per unit time, while in the present case it appears in the transition 
probability, i.e., time-integrated quantity. Anyhow, a physical meaning of 
the factor as well as of the pre-exponential factor are yet to be understood. 

For actual fusion reactions, one-dimensional treatment is obviously an 
over-simplification, which is readily understood by considering a mass-asymmetric 
entrance channel. In addition, neck formation etc. would come into play. 
It is worth to notice that even in such situations, i.e., in multi-dimensional 
problems, we can derive the same type of formula as Eq. (6), starting with an 
assumption of a quadratic potential generalized to a multi-dimension. That 
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indicates that we can define an effective one-dimensional model. Thus, the 
qualitative understandings obtained above with the schematic one-dimensional 
model are considered to be useful, with the barrier height etc. being con- 
sidered to be effective quantities. 



4. Approaching Phase ; Passing-Over Coulomb Barrier under 

Friction 

One could apply the formula obtained in the previous section to passing- 
over Coulomb barrier, approximating again the barrier as an inverted parabola. 
There, however, is another problem, as stated in section 2. In the approach- 
ing phase, dissipation of the orbital angular momentum comes into play, cou- 
pled with the radial motion. For the problem, the most simple and readily 
applicable model is the surface friction model (SFM), proposed by Gross 
and KalinowskijS] in order to explain so-called Deep-Inelastic Collisions. 

Below, we reformulate it, starting with a general framework of Ref. [20] 
which includes so-called rolling friction. Starting with intrinsic spins of the 
incident ions, L\ and L 2 , respectively, we introduce the following variables, 



L + = L x + L 2 = L - L{t) (10) 
L~ = (C 1 L 2 -C 2 L 1 )/C 1 + C 2 , 

where Lq denotes an incident orbital angular momentum and Cj, i being 1 
or 2, is an effective ion radius defined as follows, 



d = Ri(l- (b/Ri) 2 ) , (11) 



where b = lfm and Rt = 1.28 • A 1 / 3 - 0.76 + 0.8 • A i 1/3 with A; being the 
mass number of i-th ion. Then, a Langevin equation for two-body collisions 
is written as 




where fi is equal to the reduced mass of the entrance channel, and V denotes 
a sum of the Coulomb V c and the nuclear V n potentials with the rotational 
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energy given by the orbital angular momentum L. The friction tensor 0ij 
and r are given below, 



0r 



011 012 
021 022 



C r *(r)//z 
'-¥(r) • C T 

-*(r) • C ro u • ff 2 



(15) 



- 4- 

P (C\+C 2 ) 2 
1 



C1J2—C2J1 



-*(r) • C* T 



C1+C2 J1J2 



r 2 . C1J2-C2J1 " 

C1+C2 J\-J2 



Jl-h 



where Vl/(r) is a form factor specified below, and Cy and C ro u denote strengths 
for tangential and rolling frictions, respectively, In addition, a parameter g 
is introduced for describing a effective depth of the rolling friction, which is 
taken to be l.Ofm. Jj, i being 1 and 2, are the rigid moment of inertia of 
the incident ions which are assumed to be spherical. Then, the strengths r 
and 9ij are adjusted to satisfy the dissipation-fluctuation theorem with the 
friction tensor r and 0ij. The coefficient 0u is given by 0u — \l/(r) • Ct/^- 
Langevin forces are given by Wi, i being r, 1 and 2 which denote Gaussian 
random numbers and are assumed to have the following properties, 



fa) = 0, (uiWujtf)) = 2 • <% ■ 5(t - (16) 

If one wants to introduce deformations of the ions, one has to introduce 
additional degrees of freedom which describe their orientations. If we assume 
that the rolling friction is very weak compared with the others, we take C ro u 
to be zero. Then, dL — jdt = 0, and L_ = constant = L_(— oo) = 0. The 
equation for the orbital angular momentum is rewritten simply as follows, 



dL 
~dt 



= -K$/(i ■ (L - L st ) + flu • lo 1 , 



(17) 



where the effective friction and the limiting angular momentum L s t are 
given as 



Ks = C T ■ 



fir 



{Ci + C 2 f 



v J J , 

If we approximate Ri = Ci, 



'Cl Cf 

Jl + J2 



*(r) 



(18) 



- -C T *(r), 



L 



St 



(19) 
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Fig. 1. Results of SFM calculations for 48 Ca+ 244 Pu system. The top panel shows 
the sticking probability -P s ti c k as a function of incident energy relative to the 
Coulomb barrier. The middle panel shows the distribution of the radial momentum 
at the contact point, where P is given in unit of 10~ 21 sec-MeV/fm. The bottom 
panel shows the orbital angular momentum divided by the limit L st as a function 
of relative distance. 

L s t is so-called rolling limit, while one could obtain the sticking limit if one 
takes the limit that the drift part of the r.h.s. of Eq. (15) is equal to null 
vector, as discussed in Ref. [21]. Together with K r = fj,[3 r = C r ■ *(r), 
Eqs. (12), (13) and (17) just correspond to SFM. The correspondence is 
precised by giving the following relation of the friction forces, 

* (r) = (^r) 2 ' K = lc T = 0M, K? = C r = A, (20) 

where the numerical values are given in unit of 10 -23 s/MeV. The dissipation- 
fluctuation theorem is satisfied by the equations: 6> 2 = K r ■ Ta and Q\ x = 
r 2 • ■ Ta with Ta = ?a(<) being temperature of the colliding system in 
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the approaching phase. Examples of numerical solutions with the proximity 
model and with SFM are given in Ref. [22]. Generally speaking, the for- 
mer is much weaker than the latter, but for the moment we cannot draw a 
definite conclusion on which one is correct or more realistic. The former ne- 
glects frictional force stemming from strong inelastic excitations etc., while 
the latter does a rolling friction. 

We have applied SFM to superheavy systems, such as Ca + actinide 
targets. As an example, we discuss results on 48 Ca + 244 Pu system in detail. 
The top panel of Fig. 1 shows probability P s tick for the entrance system to 
reach the contact point, i.e., the relative distance being equal to a sum of 
the half density radii of the incident ions as a function of £ c . m . relative to 
the barrier height. If there is no friction, it should be always equal to 1 
above the barrier height (Below the barrier, it is equal to zero, since the 
equation is classical.). But the results are not like that. It starts with an 
extremely small value at the barrier height energy and slowly increases to 
reach 1/2 about 12 MeV above the barrier, which would explain a part of the 
extra-push. More interesting is that a distribution of the radial momentum 
calculated at the contact point is found to be approximately of a Gaussian 
one with its average value being almost exactly equal to zero as shown in the 
middle panel of Fig. 1, which indicates that the relative motion is completely 
damped at the contact point. In fact, the average orbital angular momentum 
also approaches to the dissipation limit at the contact point, which is shown 
in the bottom panel of Fig. 1. This is also the case for other actinide targets, 
say, 248 Cm and 252 Cf. 

In brief, the analyses of the approaching phase provide us with sticking 
probability -Pstick(£"c.m.) as well as with information of the amalgamated 
system, with which we can start to solve a Langevin equation for shape 
evolutions and then, can obtain formation probability -Pform- This means 
that we treat the two-body collision processes and shape evolutions of the 
united system consistently. 

5. Realistic Calculations of -Pf or m and Fusion Cross sections 

In order to describe shape evolutions starting from the pear-shape con- 
figuration of the amalgamated system to the spherical shape, at least, three 
parameters, say, distance between two mass centers R, mass asymmetry 
a, and neck parameter e in the Two-Center Parameterization |23j . With 
OBMjTH]. the friction for the neck degree of the freedom is much stronger 
than the others and thus its motion is considered to be much slower than the 
other twos. Then, we expect that two variables could describe the formation 
dynamics reasonably well, with the neck parameter e freezed. We again em- 
ploy a classical dissipation- fluctuation description, though quantum effects, 
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such as a tunneling effect, might play a significant role in passing over the 
saddle. 

A multi-dimensional Langevin equation is written as usual 24 , 

^ = (m-%- Pj (21) 

dpi dU 1 d r _ r 

~dt = ~ %~2%^ )jk-Pj-P R -Hj-(m ) jk -Pk + gij-^j, 

where qi, i being 1 or 2, specifies R, or a, and summations are implicitly 
assumed over the repeated indices. The inertia mass tensor rriij is calculated 
by Werner- Wheeler approximation 25^ and the friction tensor 7^ by OBM 
as functions of the variables R and a. The potential U is given also by the 
macroscopic LDM energy. In case of a finite angular momentum, it should 
include the rotational energy calculated with the rigid moment of inertia. 
The random force in the r.h.s. of Eq. (21) is assumed to be Gaussian and is 
expressed with a Gaussian number Tj and a strength tensor which are 
assumed to satisfy the following properties, 



(Tiit)) = (22) 
(Ti^-Tjit')) = 2-S ir S(t-t')c 

9ik9jk = lij ■ T, 

where () denotes an average over all the possible realizations. The last 
equation expresses the dissipation-fluctuation theorem. In order to obtain a 
formation probability, i.e., a probability for the system to overcome the con- 
ditional saddle point or the ridge line, we have to calculate a large number 
of trajectories. 

Examples are shown in Fig. 2, for 48 Ca- 238 U system with zero initial ra- 
dial momentum but with the temperature corresponding to the excitation 
energy 70MeV, starting at the contact configuration. It is seen that some 
trajectories go into the spherical configuration and its around, while the oth- 
ers go back to re-separation. The formers consist a formation probability, 
while the latters do quasi-fission components which are to be carefully anal- 
ysed in a future, including deformations of nascent fragments, mass drifts 
etc. before scission. 

Formation probabilities calculated with e being 0.8 are shown in the 
upper panel of Fig. 3, for 48 Ca- 238 U system. And fusion probabilities calcu- 
lated by Eq. (2), i.e., obtained by combining with the sticking probabilities 
obtained by SFM, are shown in the lower panel of Fig. 3 for the total an- 
gular momenta of the system J = and 30. Then, excitation functions of 
fusion cross sections are calculated according to the following formula, 
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Fig. 2. Examples of trajectories in two-dimensional space of the relative distance R 
divided by the radius of the sphere of the total system Rq and the mass- asymmetry 
a for 48 Ca+ 238 U system. Initial momenta are taken to be zero, but Langevin forces 
are calculated with the temperature given by the excitation energy of 70MeV. 



^fusion = S(2J + 1) • Pf us ion- (23) 

The results for 48 Ca + actinide target systems are shown in Fig. 4, compared 
with the available experimental data obtained at GSI 26 and Dubna 27 . 
It is surprising that the calculations reproduce the experimental data very 
well, not only their absolute values, but also their energy dependence, sys- 
tematically over three systems. A prediction is made for 252 C/ target case, 
which should be verified by experiment. 

6. Survival Probability 

The survival probability P S urv is a probability for the compound system 
to survive against fission decay and charged particle emission. We first 
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Fig. 3. The upper panel shows the formation probability -Pform calculated for the 
total spin J = and 30 cases of 48 Ca+ 238 U system, as a function of incident energy 
relative to the Coulomb barrier. The lower panel shows the corresponding fusion 
probability Pf us ion calculated together with P s ti c k by SFM. 



discuss the characteristic features qualitatively and then present examples of 
realistic calculations made by a new statistical code KEWPIE|17j in the next 
section. Since decay widths for the latters are small compared with those 
for the former and neutron emission, the total decay width is approximately 
given by T f + T n , and then the survival probability is given approximately 
by 

p BU rv = r n /(r / + r n )^iyr /) (24) 

where Tf and T n denote fission decay and neutron emission widths, given 
by Bohr-Wheeler [2*5] and Weisskopf 29 formulae, respectively. Of course, if 
intrinsic excitation energy E* is large enough for emissions of more than one 
neutron, the expression of Eq. (24) is repeatedly used in multiplication. In 
SHE, S> r n , so the second equation approximately holds in superheavy 
nuclei generally except cases with very large shell correction energies, which 
is easily seen by their approximate expressions, 
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Fig. 4. Calculated fusion excitation functions are shown for 48 Ca+ 238 U, + 244 Pu, 
+ 248 Cm, and + 252 Cf systems, together with the available experimental data from 
GSIH3 and Dubna|27]. 



r n ^e~ B "/ T , T f ~e- B f' T , (25) 
then the probability is given by 

Psurv =i e~( B "~ B f) /T , (26) 

where Bf and B n denote fission barrier height and neutron separation en- 
ergy, respectively. And Bf is almost equal to minus of the shell correction 
energy, because macroscopic fission barriers, i.e., LDM fission barrier Bj DM 
is very small and is nearly equal to zero for SHE, due to the fact that the 
fissility parameter Xf is close to 1. 

It is worth to consider how an excitation-energy dependence of the shell 
correction energy comes into play. As is expected, absolute values of the 
shell correction energy are reduced by excitation, so in the beginning of de- 
cay process. This is well taken into account by Ignatyuk's prescription [3U] of 
excitation-energy dependence of the level density parameter of the spherical 
shape, i.e., for neutron emission, 



a n = a n -[l + f(E*)-SE/E*}, 
f(E) = l-exp[-E*/E d ], 



(27) 
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where a n is an asymptotic level density parameter in high excitation, and 
E* intrinsic excitation energy of the compound nucleus. 5E and Ed denote 
shell correction energy of the ground state and so-called shell damping pa- 
rameter, respectively. The parameter Ed is obtained to be about 18 MeV 
by calculating excitation energy dependence of the free energy with a single 
particle model. With Eq. (27), the fission width is approximately given as 
follows, 



Tf e - B c$ /T , (28) 
B cfi = B f + f(E*) • SE. 

Then, asymptotic behaviors for E* <C Ed and E* 3> Ed become as 
follows respectively, 



B eS Bf + E* ■ 5E/E d B f , E* « E d , (29) 
B f + 5E^ Bj DM , E* > E d , 

where Bf = 5j DM — 5E denotes the fission barrier height of the ground 
state. As is seen from the above arguments, the survival probability P sur v 
is crucially determined by absolute values of the shell correction energy!! 
Remaining ambiguities are Kramers jl 11 115j and collective enhancement [To] 
factors. The former takes into account an effect of friction force acting on 
the fissioning degree of freedom, and is given by 

K f = \/x 2 + 1 - x. (30) 

This is always smaller than 1 and is approximately equal to 1/x in case 
of large x. The collective enhancement factor takes into account a differ- 
ence between collective level densities at the spherical shape and the saddle 
point shape. Since the saddle point shape of SHE is determined by shape 
dependence of the shell correction energy, no simple formula is available. It 
should be worth to notice here that so-called Strutinski correction factor [3*T] 
for Bohr- Wheeler formula fi -oj/T can be considered to be a part of the col- 
lective enhancement factor, i.e., the part from the fissioning collective degree 
of freedom, though the main part of the enhancement is expected to be that 
of the rotational degrees of freedom. 



7. Preliminary Results of Residue Cross Sections 

In order to make realistic calculation of the survival probability we have 
made a new statistical code KEWPIE (Kyoto Evaporation Width calcu- 
lation Program with time Evolution) 17_. This program treats both the 
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Fig. 5. Residue cross sections calculated by the statistical code KEWPIE arc shown 
for 48 Ca+ 208 Pb, compared with the experimental data 38 . 



production of residue as a function of the time and the final residue pro- 
duction. In the present case we will only consider the amount of nuclei 
remaining at the end of the disintegration cascade. Detailed formalism and 
the computer code will be published elsewhere. 

This program includes the main features required in the section 6. How- 
ever in this code the evaporation width of particle is calculated more ac- 
curately in the Hauser-Feshbach formalism [32]. Moreover the evaporation 
of protons, alphas and gammas are included in the program. Calculation 
of fission width is done according to Bohr- Wheeler formula with the trans- 
mission coefficient by Hill and Wheeler [6'6\ and with Strutinski correction 
factor. The fission barrier B^ DM is that of the empirical formula given for 
heavy elements by K.H. Schmidt et al. in reference |34| . 

The level density parameters a n and aj are calculated with Toke and 
Swiatecki iormulafH]^ by taking into account of shapes of the ground state 
and the saddle point. At the ground state the shape of the nucleus is assume 
to be spherical and we take into account the shell correction effect with the 
Ignatyuk prescription [3U| with Ed = 18 MeV as given in Eq. (27). At saddle 
point, deformation is evaluated by the Hasse and Myers formula [SHI and 
no shell correction effect are taken into account. 

The KEWPIE calculation has few free parameters, the scaling factor of 
the shell correction taken from M0ller et al.'s table[2] and the parameters of 
Kramer factor Kt. The latter is calculated with fiuo = 1 MeV and a friction 
factor = 5 X 10 20 sec _1 . 
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Exitation energy (MeV) 

Fig. 6. Calculated residue cross sections are shown for 48 Ca+ 244 Pu system, together 
with Dubna data 39 . Details on the parameters are given in the text. 

For 48 Ca+ 208 Pb system, fusion probabilities are calculated with the 
proximity potential 37 . because no fusion hindrance is observed there. With 
the parameters fixed, we calculate xn residue cross sectoins, whose results 
are shown in the Fig. 5. The experimental cross sections |38| are seen to be 
well reproduced, which appears to guarantee the code KEWPIE. 

For 8 Ca+ 24 Pu reaction, we use fusion probabilities calculated accord- 
ing to Eq. (2) with the realistic calculations of P{ ovm given in section 5. 
As discussed in section 6, the crucial parameter in the survival probability 
P surv is the shell correction energy. The scaling factor of 2/3 or even smaller 
has turned out to be necessary to the shell correction energies of P. M0ller 
et al. in order to be consistent with the data|39j. as is shown on Fig. 6. 
As theoretical values of the shell correction energy are very different from 
one model to another, more precise investigations are desired. We are now 
studying the reactions of 48 Ca + actinide targets, using several predictions 
of the shell correction energy, which could be very informative on the models 
of nuclear structure for heavy and superheavy nuclei. 
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